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We present recent progress in our calculation of Bk with improved staggered fermions using 
chiral extrapolations based on SU(3) staggered chiral perturbation theory. We have accumulated 
significantly higher statistics on the coarse, fine, and ultrafine MILC asqtad lattices. This leads 
to a reduction in statistical error and an improved continuum extrapolation. Our updated result 
is Bk = B/f(RGI) = 0.737 ±0.003(stat) ±0.046(sys). This is consistent with the result obtained 
using chiral extrapolations based on SU(2) staggered chiral perturbation theory, although the total 
error is somewhat larger with the SU(3) analysis. 
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1. Introduction 

Calculations of the kaon mixing parameter Bk have reached a mature stage with several re- 
sults available having all errors controlled and small. Our calculation uses improved staggered 
fermions — HYP-smeared valence on asqtad sea — and at present achieves total error of ~ 5% [|l|,0. 
Our main analysis uses chiral fitting functions derived from SU(2) staggered chiral perturbation 
theory (SChPT), and is updated in Ref. [Q]. As a check on this result, we also carry out an anal- 
ysis using fit forms from SU(3) SChPT, and in the present report we update the results from this 
analysis. In particular, we focus on the progress since our publication [|]] and last year's lattice 
proceedings 

In Table |I| we show the current set of ensembles on which we have done the SU(3) analysis, 
together with the resulting values for Bk- In the last year, we have accumulated higher statistics 
on the C2, C5, Fl, F2, and Ul ensembles. The most important improvements are those for the FI 
and Ul ensembles, since these are used for continuum extrapolation. As described in Ref. [Q], the 
9-fold increase in statistics on the Fl ensemble significantly impacts the continuum extrapolation 
in the SU(2) analysis. It forces us to fit to only the three smallest lattice spacings, excluding the 
coarse ensembles. Our main aim here is to show how the increase in statistics impacts the SU(3) 
analysis. 



a (fm) 


amt/ants 


geometry 


ID 


ens X meas 


Bk (N-BBl) 


Bk (N-BB2) 


0.12 


0.03/0.05 


20^ X 64 


CI 


564 X 1 


0.555(12) 


0.564(17) 


0.12 


0.02/0.05 


20^ X 64 


C2* 


486x9 


0.538(12) 


0.535(17) 


0.12 


0.01/0.05 


20^ X 64 


C3 


671 X 9 


0.562(6) 


0.592(14) 


0.12 


0.01/0.05 


28^ X 64 


C3-2 


275 X 8 


0.575(6) 


0.595(13) 


0.12 


0.007/0.05 


20^ X 64 


C4 


651 X 10 


0.564(5) 


0.598(13) 


0.12 


0.005/0.05 


24^ X 64 


C5* 


509x9 


0.576(5) 


0.598(12) 


0.09 


0.0062/0.031 


28^ X 96 


Fl* 


995 x9 


0.536(3) 


0.561(10) 


0.09 


0.0031/0.031 


40^ X 96 


F2* 


850 X 1 


0.539(7) 


0.544(13) 


0.06 


0.0036/0.018 


48^ X 144 


SI 


744x2 


0.535(6) 


0.560(11) 


0.06 


0.0025/0.018 


56^ X 144 


S2* 


198 X 9 


-NA- 


-NA- 


0.045 


0.0028/0.014 


64^ X 192 


Ul* 


705 X 1 


0.543(4) 


0.554(8) 



Table 1: MILC asqtad ensembles used in the calculation. Ensembles marked with a * have improved 
statistics compared to last year, while those marked with a * are new. Results for BK{fi — 2 GeV) using both 
N-BBl and N-BB2 fits are given. See the text for discussion of these fits. 



2. Fitting and Results 

In our numerical study, our lattice kaon is composed of valence (anti-)quarks with masses /Mj^ 
and rriy. On each MILC asqtad ensemble, we use 10 valence masses: 

antx, anty = ants X n/\0 with n = 1,2,3, 10 (2.1) 
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where ants is the nominal strange sea quark mass of the given lattice ensemble. Hence, we have 55 
combinations of the lattice kaon: 10 degenerate (m^ = niy) and 45 non-degenerate (nix ^ my). 

An important difference between the SU(2) and SU(3) analyses is the number of the mass 
combinations that are included. In the SU(2) analysis, we use only those in which the valence d 
quark is much lighter than the valence s quark. Specifically, we use the lightest 4 or 5 values of 
ntx and the heaviest 3 values of niy. In the SU(3) analysis, by contrast, we use all 55 combinations. 
While this makes better use of our data, it does so at considerable cost. First, quite a few of our 
mass combinations are in the regime where next-to-leading order (NLO) SU(3) ChPT is beginning 
to break down. Second, the fit forms in SU(3) SChPT contain very many fit parameters |Q] and 
have to be simplified by hand in order to be practical. In the end, as explained in Ref. [[I]], we came 
up with two different schemes for fitting, both using Bayesian constraints on parameters which 
arise due to discretization errors, but doing so in somewhat different ways. We focus here on the 
results for these two schemes, which we call "N-BBl" and "N-BB2". Both schemes are based on 
NLO SChPT with the addition of a single analytic next-to-next-to-leading order term. 

In the N-BB 1 scheme, we fit the data in two stages. First, we fit the 10 degenerate mass combi- 
nations to the functional form of Eq. (62) of Ref. [|l|] using the Bayesian method. We constrain the 
coefficient of the term arising from discretization and matching errors assuming that discretization 
errors dominate, so that C4 <x (qAqcu)^- Second, we fit all 55 combinations to the fitting functional 
form of Eq. (63) of Ref. ^ using the Bayesian method. In this second stage, we make similar as- 
sumptions concerning the size of terms arising from discretization errors, and also input the results 
from the degenerate fit (for a subset of the total set of parameters) as Bayesian constraints. 

The N-BB 2 scheme differs only in that we assume that matching errors dominate in the low- 
energy coefficients arising from discretization and matching errors, so that, for example, C4 a^, 
(where we evaluate at the scale of jU = 1/a in the MS scheme). 

In Fig. [Tj, we show how the 9-fold increase in statistics impacts the N-BBl fits on the Fl 
ensemble. As expected, errors have been reduced by a factor of ~ 3, and this holds also for the 
value for Bk that results after extrapolation to physical valence-quark masses and removal of taste- 
breaking lattice artifacts. This final value (shown as a blue octagon in the figure) has shifted up by 
about 0.2a. 

The changes in the N-BB2 fits (not shown) are different: the final value is shifted up by 1.8a 
while the error is reduced only by a factor of ~ 1 .2. 

3. Continuum extrapolation 

In Fig 1^ we present our updated results for the continuum extrapolation using both N-BB 1 
and N-BB2 fits. In the case of the N-BBl fit, it is not possible to obtain a good fit to all four values 
of a using a simple fitting functional form, e.g. ci + 020^ + ci,a'^ , with physically reasonable values 
for the coefficients. This is the same issue that arises in the SU(2) fits, and is discussed in that case 
in more detail in the companion proceedings [||]. We proceed by fitting only to the smallest three 
values of a, using either constant or linear fits. 

For the N-BB2 fit we do obtain reasonable fits to all four values of a, but, in order to compare 
with the N-BBl fits, we also use only the smallest three values of a. 
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Figure 1: Bk{1 /a) versus Xp (squared mass of pion composed of valence x and x) for the Fl ensemble, 
using N-BBl fits. The left panel uses 1 measurement/configuration, while the right panel uses 9 measure- 
ments/configuration. Red diamonds show the data, while the blue octagon shows the result obtained after 
extrapolation to physical quark masses with all taste-breaking lattice artifacts removed. 
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Figure 2: B/f (NDR, ^ = 2 GeV) as a function of a^ (in fm x 100) for the N-BBl fit (left) and the N-BB2 fit 
(right), showing constant and hnear extrapolations to the smallest three values of a. 



In Table ^, we summarize the results of these continuum extrapolations. We use the constant 
fit to the N-BBl data for our central value and the difference between this and the result from 
the corresponding fit to the N-BB2 data as an estimate of the fitting systematic. Note that the 
difference between N-BBl and N-BB2 fits is much larger than the difference between constant and 
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fit type 


constant fit 


linear fit 


N-BBl 
N-BB2 


0.5383(23) 
0.5573(53) 


0.5344(53) 
0.5526(106) 



Table 2: Results for = 2 GeV to data from N-BBl and N-BB2 chiral fits using both constant and 

linear continuum extrapolations (in a^). 



linear extrapolations. This reflects the uncertainty in SU(3) fitting caused by the large number of 
parameters related to discretization and matching errors. 

In Fig. ^ we compare the results from the N-BBl fits to those using our preferred SU(2) fitting 
approach. There is reasonable consistency point by point. 
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Figure 3: Bj? (NDR,^ = 2 GeV) as a function of a^ (in fm x 100) from the SU(3) analysis (N-BBl fit) in 
the left panel and from the SU(2) analysis in the right panel. Constant fits to the three smallest values of 
a = are shown. 



4. Error Budget and Conclusions 

In Table H we list various sources of error in the SU(3)-based calculation of Bk- Many of the 
smaller errors are unchanged from Ref. and we refer to that paper for an explanation of how 
they are estimated. 

The major changes to the errors are as follows. The statistical error has been reduced from 
1.4% in Ref. ^ to 0.4%, due to the use of more measurements. The "matching" error — due to 
our use of one-loop matching between lattice and continuum operators — has been reduced from 
5.5% in Ref. [|l|] to 4.4%. This is simply due to the addition of the smallest lattice spacing, for our 
estimate of the percentage error is (/x = l/amin)^- The discretization error has also been reduced 
(from 2.2%), since our smallest value of a is closer to the continuum. Finally, the "fitting (2)" 
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cause 


error (%) 


memo 


status 


statistics 


0.43 


N-BBl fit 


update 


matching factor 


4.4 


AB?^ (Ul) 


update 


discretization 


0.95 


diff. of const and linear extrap 


update 


fitting (1) 


0.36 


diff. of N-BB 1 and N-B 1 (C3) 


m 


fitting (2) 


3.5 


diff. of N-BB 1 and N-BB2 at a = 


update 


ami extrap 


1.0 


diff. of (C3) and linear extrap 


m 


arris extrap 


0.5 


constant vs. linear extrap 


m 


finite volume 


2.3 


diff. of 20^ (C3) and 28^^ (C3-2) 


® 


n 


0.12 


error propagation from ri 


t 



Table 3: Error budget for Bk obtained using SU(3) SChPT fitting. 



error — that due to the uncertainty in the size of SChPT coefficients introduced by discretization 
and matching errors — has changed from 5.3% to 3.5%. Nevertheless, this error remains large, and 
along with the matching error, dominates the total error. This large fitting systematic is a reflection 
of the difficulties in using SU(3) SChPT and is reason why we think that the SU(2) analysis is more 
reliable. 

Our present value from the SU(3) analysis is 

B/f (NDR, n=2 GeV) = 0.5383 ± 0.0023 ± 0.0337 , SU(3) fit 
Bk = Bk{RGI) =0.7371 ±0.0032 ± 0.0461, SU(3) fit 

where the first error is statistical and the second systematic. The total error is 6.3%. This should be 
compared to our updated SU(2) result 

Sjf = B^(RGI) =0.725 ±0.004(stat)±0.038(sys) SU(2) fit. (4.2) 

We see that, although the SU(3) result has a smaller statistical error, it has a significantly larger 
systematic error. The most important observation, however, is that the two results are consistent. 
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